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1. PURPOSE 

1.1 To compute in floating point the modified Bessel function Io (x) 

2. RESTRICTIONS 

2.1 Require AN-044 exponertial subroutine, 

3. METHOD 

3.1 Truncated power series or exponential approximation, depending on size of 
argument. 

3.1.1 X<7: Io<*l + (^> 2 + jr (f) 2(k+1) 

k = 1 G k 

3.1.2 X2.7: Io^ <£ J~l + Jx + Jr*2 + §| , 3 ) 

3.2 Accuracy within £ 0.01 % 

3.3 Reference 

H, 3. Dwight Tables of Integrals and Other Mathematical Data 
MacMillan Co., 1947. Equation 813.1 and 814,1. 

4. USAGE 

4.1 The magnitude of the argument must be in the A and R register as a positive 
normalized floating point number. The result is left in the A and R register 
upon exiting. 

4.2 Calling sequences 

4.2.1 If AN 044 is located at (Io + 140.0) 

TRA Io 
Return 

4.2.2 If AN 044 is elsewhere in memory 

TRA (Io + 132.0 octal) 
PZE EX£ 



e<_ + 1 Return 
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EXP = first location of AN-044 

Io = first location of Io routine. 

After you have used 4.2.2 on a run the program can now be 
entered by 4.2.1. 

4.3 Uses locations 0000 - 0136 relocatable. No temporary storage. 

4.4 Brror Codes 

If the limits of AN-044 are exceeded there will be a halt and transfer 

to (Io + 0043.1). If you press start you will return to the main program. 

4.5 Both L and V loops are used by this routine. 
5, Coding Information 

5.1 Constants (floating point) 



2 tr 



0100 
0102 
0104 
0106 
0110 
0112 
0114 
0116 
0120 
0122 
0124 
0126 
0130 



7.0 

6.283185308 

1.0 

4.0 

36.0 

576.0 

14400.0 

518400.0 

23401600.0 

2.0. 

37.5 

8.0 

4,5 
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0000. +15.0070.0 — SAX 



5.0070. 
IT0OO7; 



+01.0007.0 ADD 



0001 .0 +H2.OOH3.I STA 
+ 1|^OOOOJO— JCAR. 



0002.0 +60,0071 .0 — SIQ 



+30.0070.0 FCA 



0003.0 +06.0100.0 FSB 
+ 52 r 0QWl,O.--TRP- 



000>l,0 +6U.001Q.0 — CXL 



+66.0020.0 CTV 



0005.0 +30.0070.0 FCA 

+05.^01 -22.a.__FDy. 



0006 ,o +35.777^ .n — ejsi 



+57.776O.O TRA 



0007.0 +00.0000.0 

— **QQ»Q0QQ.4— «=■- 



.-ooi-a»o-4O7*7HP0i*a-XMP. 

+35-777^.0 FST 



0011.0 +O7.777H.O FMP 
*35> 7776.0 FST 



._GG4£-^-+a5_cuii6.o___EDy_ 

+04.7772.0 FAD 



0013.0 +35.7772.0 FST 
+30.777^.0 FCA 



._0CUiU0-+-Q7.77-76.0___M > _ 
+35.7776.0 FST 



0015.0 +05.0110»0 FDV 
+ pfi , 7772 , FAR 



-9€H-6 i r0-+35w77 L 72»O— £ST- 
+30.777^.0 FCA 



0017.0 +07.7776.0 FMP 
I - 35.7776,Q FST 



-O02O.-O -+0R- r 04 A -2*0— FDV . 
+04.7772.0 FAD 



0021 .0 +6H. 0030.0 CTL 
+5.7^7750^0— TRA - 



-0022*0 _+oa^oooo^o— , ™= 
-00.0000.0 — 



0023 .0 +00.0000.0 Page k 

-ofr.-eoQO-.-o--—- 



-GG2+-.-G -+ oa.-oooa^o - - — 

-00.0000.0 



0025.0 +00.0000.0 

^oaux>o£U}-=r----_ 



-00.0000.0 — 



0027.0 +00.0000.0 — 
OGs-0000,0- 



0030.0 + yj. 7770.0 — fst- 

+30.777^.0 FCA 



0031.0 +07.7776.0 FMP 
+-35.-r??6;0— FST- 



0032.0 +05.0H 4.0 — FW- 
+04. 7772.0 FAD 

0033.0 +35.777,2.0 FST 



003*1. +07 "7776*0 FMP 



+35.7776.0 FST 



00^5.0 +O5.OI 16.0 FDV 
_ — +Q*U777-2*Q- - -F-AD. 



0036,0 +35- 777 P. O FST 
+16.777^.0 FCA 



0037". ~+Gk ."66¥6."6~ " "cf l 
+5J.77-&.0._.IRA. 



0040. +07.7776.0 FMP 
+^.-7776;0---FST- 



U04l .U +US.0I20.0 — FW 



>4. 



+04.7772.0 FAD 



0042.0 +04.0104.0 FAD 

+t)if;777 i f;0r-rAcr 



00+3.0 +4U.UUUU.0 — AttS" 
+57.2001 .0 TRA 



0044.0 +30.0070.0 FCA 

+57;0T+tr,xr-TRA 



0045.0 +00.///4.0 CLA 
+77.0043.1 HTR 



00+6.0 +64.0050.0 CTL 
+6^70060.-0--- CTV 



00^7.0+35-7776.0 FST Page 5 

+ l 57. 7760.0 TRA 



00*50.0 +^0.0102.0 FCA 
— - +O7.OO7O.O — FMP- 



t)0^ T -."O -+^5V777^."0 " " FST " 
+44.7774.0 FSQ 



0052.0 +35-7774.0 FST 

+^0. 7 7 7 6.0 — fc*- 



"0053';o-+tK:777^".^"-i ? ^ 

+35.7776.O FST 



005^.0 +30.01 04.0 FCA 

+05.0126.0 — rw- 

"0055".~0- +-05: 0070-.0-- ~rm 

+35.777^.0 FST 



OO56.O +07.7774.0 FMP 
+=i5.0Q 7 6.Q FST 



005T.~0" +O7.~0130; tr - -FMP- • 
+04.7774.0 FAD 



0060.0 +35.0074.0 FST 



~oo6r.-o-T07;0076iO--FMP- 

+07.0124.0 FMP 



0062.0 +04.0074.0 FAD 
+04^01-04^0— FAD- 



-0063^0-^07^-7776-.^— fmp- 

+57.0043.1 TRA 



oo64.o +00.0000.0 — 
---oe^oooe^€>— — 



-0065^0- +--00i0000^0— -—■ 
-00.0000.0 



0066.0 +00.0000.0 — 
--oo;0000;0— -— - 



-oo67;0-^oo;oooo-.o— -*-- 

-00.0000.0 — 



0070.0 +40.0000.0 ARS 



00 7 1.0+00.0000.0 — 
-00.0002.0 



0072.0 +00.0000.0 

-oa.-ooooro— «•** 



0073-0 +00,0000.0 Page Six 

-00.0000.0 — 



007*+. 4-^2.2000.0 ST A 

~--ooygo©oyO--^----=- 



00 7* 5.0 -00.0000.0 — — 
-00.0002.1 ~ 



0076.0 440.0000.0 ARS 

-oo-.-oeoo-.-o--^---^- 



uu/7.0 -00. oooo. — =■ 

-00.000*5.1 — 



0100.0 4-70.0000.0 

-uu.uuuu.o — =■ 



"otot;o"4-uo:ouoo".o— --- 

-00*0001.1 



0102.0 4-62.2077.0 

+52>AO.U — TRF 



"OTOT.'O" +do;oooo;o— --- 
-00.0001.1 — 



010*1.0 +ho. 0000.0 ARS 
-00.0000.0 — 



-otO5-.-0-4-^3Oi-0OOOi-e — * — 

-00,0000.1 



0106.0 440.0000.0 ARS 

-00.0000.0 — — 



i3f&7vQ -4-OQ-.-0OO0 ^9- ~ 
-00.0001 .1 



0110.0 4-^.0000.0 FSQ 
=.-00^0000^0--™-- 



"OT1T.-0-*OOjOOOOj©--- -** - 

-00.0003.0 — 



01 1 2.0 44*+. 0000.0 FSQ 

-t?o:oooo;o— -— - 



"on3";o"+oo;oooo;o""--- 
-00.0005.0 — 



011*1.0 4-70.2000.0 

-ocr."oooa."0~~-~-- 



or t 5.13 ~4-oo.~000a.-0- - --- 
-00.0007.0 — 



0116.0 4-77*2200.0 HTR 

-06.-0000.-0- - ----- 



vr t 770 --foa.~oooo.-e- - «-*■ 

-00.0011.1 — 
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"ar?D".tJ-+ftr.-3Hfe-.-o— ST6 

-00.0000.0 

0121 .0 +00.0000.0 

_ 00v00t ^ 1 ___ ;wl , 



~or 22vo -+ Htr.-oooo-.-o — t^rs • 
-00.0000.0 — 



0123.0 +00.0000.0 — 
-t)o^ooot-^o— -=-- 



n0T2ir.-0-+^5-.-HoO0^O— FNH 
-00.0000.0 — 



0125.0 +00.0000.0 — - 

--- ----- ~ -00.-0003 ;0— — «- 



-orz6.-o-+te.-oooo;0— ars- 

-00.0000.0 -— 



0127.0 +00.0000.0 — 
- - "- — — -OQ-.~00O2;O-~- *— - 



1 30.0 +kh. 0000.0 — Fee- 

-00.0000.0 — 



om .0 +00.0000.0 — 
— -— — --ooiooai-,1— ' — ■- 



0132.0 +1^.0070*0 — 5** 
+^2.013^.1 STA 



0133.0 +01 .on6.o ADD 

-— +-^2^00*3^— STA- 



013*1.0 +ho. 0000.0 — aks- 

+00.0000.0 CLA 



0135.0 +H2.0CM.1 STA 
— 1 +5-7^0001-^1— TRA- 



0136.0 +00.0000.0 — ==■ 

-00.0001.0 



